二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率,对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序,该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想,验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合,说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动。为了提高计算规模及计算速度,作者对程序进行了并行化并针对微机机群系统进行了并行优化。实际计算表明本程序具有较高的并行效率。      

爆炸波与物体干扰流场的数值模拟

本文构造一类新型的计算效率较高的ＥＮＯ差分格式。利用该格式非定常Ｅｕｌｅｒ方程出发,对爆炸波扫过物体时形成的二维和三维干扰流场进行了数值模拟,计算结果与实验数据符合较好。     

被动式引射器内两相流数值模拟

把气相二维欧拉方程的隐式矢通量分裂近似因子分解法与固相粒子的特征线解法进行耦合,对被动式引射器内两相流场进行数值模拟,确定产生激波的喷管出口压强,分析不同引射器结构对引射效果的影响,并获得了一种最佳结构。研究结果可以为被动式引射器的设计与使用提供依据。     

两种多重网格法求解径流叶轮机械可压流动

为了加快数值求解叶轮机械内部流场的收敛速度,并改善计算结果,在数值求解径流叶轮机械可压流动的计算中,采用了由微元、二重和三重网格逐级转入和同时转入的两种多重网格法。对德国宇航局Ｋｒａｉｎ Ｈ设计的离心压气机内部流场的数值计算表明,两种多重网格法的计算结果差别不大,尽管在收敛速度方面存在差异,但均与测量数据吻合比较好。     

正反顺序加载下安装边螺栓预紧力分散度的数值模拟与试验验证

在航空发动机机匣安装边装配过程中,不同位置螺栓之间会产生弹性交互作用,不同加载方案会使装配后的螺栓预紧力存在一定的分散度。针对这一问题,建立了3层带止口机匣安装边3维接触有限元模型,研究不同加载方案下螺栓预紧力的施加方法。基于有限元模型数值模拟了3层机匣安装边分别在顺序加载方案、交叉加载方案和正反顺序2轮加载方案下螺栓预紧力的分散度。搭建了相应的试验装置,分别对3种加载方案的预紧力分散度进行试验验证。结果表明：采用正反顺序2轮加载后的螺栓预紧力分散度最小,在3层止口机匣安装边上采用正反顺序2轮加载方案后的螺栓预紧力分散度较小,在有限元仿真中螺栓预紧力分散度为2.43%,在试验中螺栓预紧力分散度为4.98%。采用此方案加载后,止口机匣安装边能达到更好的密封效果和较好的刚度对称性。     

高速旋转冲压增程弹用进气道复杂流场数值模拟

采用块结构网格与二阶精度流场分区求解技术,对高速旋转、含侧向支柱冲压增程弹丸进气道内外复杂流场进行了数值模拟,得到了高速旋转工况下对应于不同来流马赫数和攻角,以及临界工况时超音速进气道内外流场复杂的波系结构。在旋转工况下,进气道流场结构与不旋转时相似,但旋转速度以及攻角的存在对冲压发动机进气道的总体性能产生了负面影响,进气道总压恢复系数和动能系数均有所降低,而流场畸变指数则显著增大;冲压发动机进气道在较低马赫数工作时的综合性能优于在较高马赫数工作时的综合性能。     

采用改进的CE/SE方法模拟方管中氢氧爆轰波的稳定传播结构

采用改进的时空守恒元/解元（CE/SE）格式和新型二步化学反应模型，对小扰动情况下方管通道中的氢氧气相爆轰传播过程进行了数值模拟，对不同截面尺寸和两种扰动模式下的爆轰波三波线结构进行了详细计算和分析。数值模拟花费了相对较少的计算代价，清晰地得到了方形截面通道中3种相对稳定的爆轰波传播结构，即直角模式、对角模式和螺旋模式。随着截面尺寸的减小，爆轰波胞格图样与波阵面结构随之产生适应性的变形，乃至出现胞格数目的减少和部分三波线的消失；最终，直角模式和对角模式都会转化为顺时针或逆时针的螺旋模式；临界尺寸下，得到了转换过程中新的过渡结构和压力脉冲振幅的变化，并以此初步讨论了直角和对角传播模式的稳定性。     

基于带波动算子非线性Schr(o)dinger方程数值分析的守恒差分算法

对一类带波动算子的非线性Schr(o)dinger方程进行了数值分析,提出了一个含参数的二阶守恒差分格式,根据参数选取的差异,该格式既可隐式计算也可显式计算.对初值条件进行了中心差分离散,使其具有二阶精度,从而与守恒格式的精度一致.利用矩阵理论证明了差分解的存在惟一性,并利用一个重要的不等式在先验估计的基础上,运用能量估计的方法证明了该格式按无穷范数以二阶精度收敛到真实解.数值实验表明该格式具有较高的计算效率.     

一种改进的MUSCL格式

给出了一种改进的二阶精度的ＭＵＳＣＬ格式，并用来求解Ｅｕｌｅｒ方程，改进后的ＭＵＳＣＬ格式在线化插值步中，用Ｖａｎ Ｌｅｅｒ的插值斜率对特片变量进行插值，并对与熵波、剪切波对应的特片变量使用了小耗散的插值斜率，在迎风步采用Ｒｏｅ的近似黎曼解求得迎风通量。比较一维和二维的计算结果均可以看出，改进后的ＭＵＳＣＬ格式在消除数值振荡的同时，明显提高了流场间断的分辨率。